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ABSTRACT 

We fully characterize the Mg II absorbing circumgalactic medium (CGM) using the "Mg II Absorber-Galaxy 
Catalog" (MAGllCAT), a compilation comprising 169 intermediate redshift (0.1 < z < 1.1) galaxies within 
200 kpc (projected) of background quasars. We examine the rest-frame equivalent widths, W r (2196), and 
covering fractions, f c , down to W,(2796) ~ 0.003 A, against galaxy redshift, B- and Zf-band luminosity, B-K 
color, and impact parameter, D. We also present B- and K-bwd luminosity functions for different W r (2796) and 
redshift regimes. The radial extent of the Mgll absorbing CGM, R(L) = R^L/L*) 13 , exhibits little dependency 
on W r (2796) threshold in the B-band, but becomes less sensitive to luminosity in the A'-band. In both the B- and 
/iT-bands, R{L) has a much steeper luminosity dependence for bluer galaxies than for redder galaxies, and for 
lower redshift galaxies compared to higher redshift galaxies. The covering fractions exhibit clear differential 
dependencies on W r (2796) threshold, galaxy luminosity, and redshift as a function of impact parameter. Lower 
luminosity galaxies exhibit no Mgll absorption beyond 100 kpc, whereas higher luminosity galaxies have 
0.1 < / c < 0.3 beyond 100 kpc, with f c decreasing with increasing W r (2796) threshold. Redshift evolution 
is found in the B-band luminosity function for galaxies associated with W, (2796) < 0.3 A absorption. Many 
of our findings are in conflict with previously published works. We address these conflicts and discuss the 
implications of our results for the low-ionization, intermediate redshift CGM. 
Subject headings: galaxies: halos — quasars: absorption lines 



1. INTRODUCTION 

Understanding the formation and evolution of galaxies 
is one of the foremost problems facing extragalactic re- 
search. Over the last decade, theoretical activity has been 
focused on studying how galaxies form and evolve in the 
context of the response of baryonic gas to dark matter ha- 
los of various masses and to ph ysical processes such as 
stellar and AGN feedback (e.g.. iBirnboim & Dekell 120031: 
Mailer & Bullocld2004llKeres et : al.ll2005blDekel & Bir nboim 
2006tlBirnboim et. al.l2007llOcvirk. Pichon. & Tevssierl2008l: 
Keres et al.1 120091: iQ ppenheimer et al. 2010; Stewa rtet al.1 
201 Ulvan de Voort et al.ll201 Ulvan de Voort & Schavd201 if . 
These theoretical works have established that the circum- 
galactic medium (CGM) is a dynamic environment compris- 
ing inflowing accretion, outflowing winds, and gas that recy- 
cles between the halo and galaxy. 

What has become clear is that the CGM is intimately linked 
to galaxy morphology, stellar populations and kinematics, 
and chemical evolution. The detailed physics governing the 
heating and cooling of baryonic gas regulates the formation 
and galactic-scale dynamical motions of stars, which in turn 
feedback and gover n the physics of the gas in th e CGM (e.g., 
Dekel &Silkl 119861 ICeverino & Kkvpir] 120091) . Since this 
paradigm of galaxy evolution suggests a strong CGM-galaxy 
connection, many investigators have taken a phenomenologi- 
cal approach to developing models that predict the observable 
geometric distribution a nd kinematics of baryonic gas in 
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As such, the CGM is an astrophysical environment which, 
at any point in the evolution of a galaxy, yields clues to both 
its historical development and future evolution. Therefore, 
observations of the CGM around individual galaxies promise 
to provide highly sought constraints on the physics governing 
the global properties of galaxies. In practice, the CGM is al- 
most exclusively accessible via absorption lines present in the 
spectra of background objects (traditionally using the tech- 
nique of quasar absorption lines). Studying the CGM-galaxy 
connection is desirable over all redshifts, however, the vast 
maj ority of work to date has been focused on z <: 1 (however, 
cf.. lAdelberger et al.1 120051: ISimcoe et all 120061: ISteidel et all 
2010), where galaxy photometric and spectral properties can 
be measured with relatively high accuracy and modest invest- 
ment in telescope time. 

Campaigns to study the CG M in absorption at z < 1 have 
targeted Ovi ab sorption (e.g.. | Tumlinson et al.1 1201 ll) . Civ 
absorption (e.g., IChen et al.1 12001 a|). and neut ra l hydrogen 
Lya absorption (e.g., Lanzetta et al.1 119951: IChen et all 
2001b). However, by far the vast majority of surveys 
have focused on Mgll a b sorpti on (iBergeron & Boisse 
199U iBechtold & E llingson l!992t 
LeBrunetalJ 1 19931 ISteidel. Dickinsc 

Churchill et all 119961: ISteidel et al.1 
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The Mgll AA2796,2803 doublet is especially well suited at 
intermediate redshifts since the absorption is observable in the 
optical. Furthermore, Mg II absorption is ideal for tracing low 
ionization, metal-enriched CGM gas in that it arises in a wide 
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range of H I environments, from logjV(H l) ~ 16.5 to greater 
than 21.5 (Bergeron & Stasiriska jl986t ISteidel & S argent 
1992t IChurchill et all 1 19991 l2000ai i Rao & Turnshekl 120001: 
Rigbv, Charlton. & Churchilll 12003 . It is also well estab- 



lished that Mgll a bsorption probes key CGM structures 
such as outflow gas (iTremonti et al. 2007; W einer et al.ll2009l : 
iMartin & Bouchel l2009t iRubin et all 1201 Ol) and inflow gas 
dRibaudo et alj 1201 U iRubin et alj 1201 ll iThom etail 1201 li 
Kacprzak et alJl2012l) associated with galaxies. 

Through the efforts of the aforementioned studies, the 
behavior of Mgll CGM absorption in relation to various 
galaxy properties has been incrementally characterized over 
the last two decades. Examples of some of the findings 
include relationships betwe en the strength of Mgll absorp- 
tion an d galaxy luminosity (Kacprzak et al. 2008: IChen et all 
2010al), galaxy mass jBouche et al. 2006; Gauthie ret alJ 
200^IChurchill et alJ 12012b!) star formation or specifi c 



star formatio n rate dChen et alJ | 2010bt iMenard et alj|201 ll) . 
galaxy color (Zibetti et al. 2007; Bo rdoloi et al 1201 lb . galaxy 



( Bordoloi et all 120111; IBouche et all 120111; 


Kacprzak et al. 


201 lbtlKacprzak. Churchill. & Nielsenll2012l; 
2012afl. 
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Generally speaking, these investigations have yielded re- 
sults that appear to be converging on an observational por- 
trait of the Mg II absorbing CGM that is more or less consis- 
tent with the broad scenario forwarded by theory and sim- 
ulations. For example, evidence is mounting that accre- 
tion may preferen t ially reside in a coplanar geometry (e.g. , 
Steidel et al.lfl997l: iKacnrzak et al.ll2010l IStewart et alJl201ll 



Kacpr zak et alj 201 lbl) , whereas winds may outflow along 
the galaxy minor axis ([Bordoloi et al J 1201 it IBouche et al. 
201 It IChurchill et"ai1l2012atlKacprzak. Churchill. & N ielsen 



20121) . However, due to the complex ionization structure, non- 
uniform metal enrichment, and dynamical processes of the 
CGM, and due to the wide range of galaxy properties, such 
as luminosity, color, mass, and morphology, the data exhibit 
substantial scatter. As such, while some works have found sta- 
tistically significant correlations between two quantities or be- 
tween combined/scaled quantities, many works have reported 
connections between galaxy and Mg II absorption CGM prop- 
erties either based on general trends (i.e., correlations that are 
not statistically significant above the 3 a level) or on slight 
reductions in the scatter of certain relationships by combin- 
ing/scaling one or more measured quantity. 

Motivated by the potential that combining the data from 
our work and other surveys may further illuminate the 
CGM-galaxy connection, we have endeavored to assemble 
a database of the extant works focused on Mgll absorption 
from the CGM of intermediate redshift galaxies. We also aim 
to provide a large uniform data set based upon a single cos- 
mological parameter set and to consolidate the measurements 
of given galaxy-absorber pairs duplicated in various works. 
Such a compilation holds the promise of yielding higher sta- 
tistical significance in the already published results, and of 
providing greater leverage for exploring the dependence of 
CGM Mg II absorption with, for example, redshift, galaxy lu- 
minosity/color, or impact parameter as a function of equiva- 
lent width threshold. 

The data and full details of the "Mgll Absorber-Galaxy 
Catalog" (MAGllCAT) we have compiled, which we briefly 
describe in § wil l be presented in Paper I of this series 
dNielsen et aUl2013l) . In this paper, we utilize the data to 
address several long-standing questions with regard to the 



Mg II CGM-galaxy connection. In § [2] we briefly describe 
the galaxy sample, which has uniform photometric proper- 
ties and impact parameters for a ACDM cosmology {Ho = 70 
km s" 1 Mpc" 1 , VIm = 0.3 and £1a = 0.7). In §[3] we present an 
examination of (1) galaxy color and absolute magnitude evo- 
lution, (2) the dependence of W, (2796) on galaxy color, (3) the 
dependence of W r (2796) on impact parameter, (4) covering 
fraction as a function of impact parameter and W r (2796) for 
luminosity, color, and redshift subsamples, (5) the luminosity 
scaling of "halo absorption radius" and the behavior of cov- 
ering fraction with W,-(2796), galaxy color, and redshift, and 
(6) B- and /f-band luminosity functions. In § [H we discuss 
the implications of our results and give concluding remarks in 

§13 

2. GALAXY SAMPLE AND SUBSAMPLES 

From an exhaustive literature search, we compiled a cat- 
alog of galaxies with spectroscopic redshifts 0.1 < z < 1 
within a projected distance of D < 200 kpc from a back- 
ground quasar, regardless of known Mg II absorption. We did 
not include galaxies for which W, (2796) measurements could 
not be obtained. We also distinguish between isolated and 
group galaxies, where group galaxies have a nearest neigh- 
bor within 100 kpc and have a velocity separation no greater 
than 500 km s" 1 . The total catalog consists of 205 galaxies. 
Full details of all galaxies and Mg II absorption properties will 
be presented in Paper I. Here we briefly describe the general 
properties of MAGllCAT. 

The sample presented here comprises 169 isolated galaxies 
toward 117 sightlines, covering a redshift range 0.072 < z ga i < 
1.120, with a median of (z) = 0.346. The impact parameter 
range is 5.5 < D < 193 kpc, with median (D) = 48 kpc. 

For each galaxy, we determined rest-frame AB ab- 
solute magnitudes, Mg and Mk, by computing k- 
corre ctions appropriate for the ob served magnitudes 
(e.g., Kim, Goobar, & Perlmutter 1996) using the extended 
IColeman et alJ (1 19801) spectral energy distribution (SED) 
templates from Bolzon ella et alj d2000t) . We selected a 
galaxy SED by comparing the observed galaxy colors to 
the redshifted SEDs. Luminosities, Lg/L B and LxjL* K , 
were obtained u sing a linear fit to Mg with redshift from 
iFaber et ail d2007l) IB-band) and using M£(z) as expressed in 
Eq. 2 from lCirasuolo et ail (1201 Oh (K-band). Absolute fi-band 
magnitudes range from -16.1 > Mg > -22.8, corresponding 
to luminosities of 0.017 < Lg/L* B < 7.44, with median 
Lg/Lg = 0.526. Absolute /T-band magnitudes range from 
-17.0 > M K > -24.3, corresponding to 0.006 < L K /L* K < 4.9. 
Rest-frame B-K colors range from 0.04 < B-K < 4.09, 
with median B-K = 1.39. We obtain B-K colors for all but 
21 of the galaxies. 

MAGllCAT was bifurcated into several subsamples for 
analysis. When dividing the galaxy-absorber pairs into sub- 
samples, we either take a sample driven approach by splitting 
the sample at the median value of a given observed quantity, 
or dividing the sample based upon historical precedent, such 
as W,(2796) cuts. The full sample was bifurcated based on 
the median galaxy redshift, L B /L B , L K /L* K , or B-K. Table Q] 
presents the characteristics of the full sample and each sub- 
sample, including the number of galaxies, the median value 
by which the catalog was cut, and the minimum, maximum, 
and mean values of galaxy redshift, B- and K-bwd luminos- 
ity, and B-K color. The subsample names listed will be used 
throughout this paper. 
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TABLE 1 
Subsample Characteristics 



Subsample # Gals Cut a 



Min 



-Zgal - 
Max 



Mean Min 



Max 



Mean Min" 



■ l k/L* k 

Max b Mean b 



B-K 

Min b Max b Mean b 



AU 


169 




0.072 


1.120 


0.415 


0.017 


7.441 


0.780 


0.006 


4.901 


0.653 


0.038 


4.090 


1.471 


Low z 


84 


0.346 


0.072 


0.340 


0.222 


0.017 


7.441 


0.740 


0.006 


4.901 


0.621 


0.430 


3.042 


1.571 


Highz 


85 


0.346 


0.346 


1.120 


0.606 


0.071 


4.801 


0.820 


0.016 


3.152 


0.685 


0.038 


4.090 


1.370 


Low Lg/Lg 


84 


0.526 


0.072 


0.941 


0.395 


0.017 


0.523 


0.276 


0.006 


1.680 


0.266 


0.038 


4.090 


1.446 


High L B /L* 


85 


0.526 


0.096 


1.120 


0.436 


0.526 


7.441 


1.278 


0.214 


4.901 


1.020 


0.483 


2.537 


1.494 


Low Lk/L* k 


74 


0.414 


0.110 


0.941 


0.386 


0.019 


1.214 


0.324 


0.006 


0.406 


0.193 


0.038 


3.037 


1.230 


High L K /L* K 


74 


0.414 


0.096 


1.017 


0.419 


0.189 


7.441 


1.231 


0.423 


4.901 


1.114 


0.483 


4.090 


1.711 


Blue 


74 


1.391 


0.098 


1.017 


0.443 


0.019 


7.441 


0.844 


0.006 


2.985 


0.443 


0.038 


1.384 


0.996 


Red 


74 


1.391 


0.096 


0.852 


0.362 


0.035 


3.759 


0.711 


0.026 


4.901 


0.864 


1.399 


4.090 


1.945 



a The median value by which each subsample was bifurcated. The cut is inclusive to the "high" and red bins. 
b Only including galaxies with a value for Mk ■ 



3. RESULTS 

In this section, we report on the luminous properties of the 
galaxies and directly compare them to the properties of the 
Mgll absorbing CGM. The analysis presented here is based 
upon direct observables. We do not scale or combine mea- 
sured quantities. We also do not scale or fit the data to models 
of the CGM. Our aim is to characterize the Mg II CGM-galaxy 
connection directly with no assumptions and to provide infor- 
mation that can be directly interpreted in terms of primary 
observables. We examine only bivariate relationships, reserv- 
ing further analysis, such as multivariate techniques for future 
work to appear in later papers of this series. 

3.1. Galaxy Redshifts, Colors, and Magnitudes 

Since the CGM as probed by Mgll may depend upon 
galaxy stellar populations, and these populations are known 
to evolve, we examined whether the galaxy rest-frame col- 
ors and magnitudes evolve with redshift for the sample. A 
Kendall-r rank correlation test on B-K versus z ga i shows that 
the null hypothesis of no correlation can be ruled out to a con- 
fidence level (CL) no better than 2.3 a 4 . A rank correlation 
test yielded a 4.4 a significance that M B correlates with z ga i, 
w hereas Mk i s cons istent with no evolution (2.1 a). 

iFaber et ail (120071) found that Mg brightens with increasing 
redshift from the DEEP2+COMBO-17 surveys. The average 
Mb of the sample presented here increases with re dshift more 
rapidly thanM|. Similarly. ICirasuolo etai1(l2010l) found that 
M K tends to brighten with increasing redshift from UKIDSS, 
whereas the average Mk of the sample presented here does not 
appear to evolve. As such, we may be seeing that galaxies se- 
lected by Mg II absorption are characterized by bluer colors at 
higher redshift. However, selection effects may be important. 
Virtually all the galaxies were selected or discovered using the 
red band in the observer frame, which probes further toward 
the B-band in the galaxy rest frame with increasing redshift. 
However, this cannot explain the lack of evolution in the K- 
band. 

Using statistical methods in which i ndividual galaxies w ere 
not directly identified with absorbers, Zibett retaTl (120071) re- 
ported no color evolution for 0.4 < z < 1 using observed g, 
r, i, and z band magnitudes for W r (2796) > 0.8 A. A direct 
comparison is difficult because we use optical B-band and 
2.2/im infrared absolute magnitudes to determine rest-frame 

4 Running the test on B—R versus z ga i gives a less significant result (1.9 ex) 



colors and we study systems with much smaller W r (2796) than 
represented by their sample. However, bifurcating the sam- 
ple presented here at W f (2796) = 0.8 A, we obtain 0.8 a for 
W r (2796) > 0.8 A. For W r (2796) < 0.8 A we ob tain 2.5 a in- 
dicatin g that the systems below the sensitivity of Zibet ti etldl 
(2007) dominate a possible color evolution. 

3.2. W r (2796), Galaxy Luminosity, and Colors 

To see whether W r (2796) exhibits some dependency on 
galaxy luminosity, we performed a non-parametric Kendall's 
r rank correlation test that allows for upper limits on ei- 
ther the dependent or independent variable for bivariate data 
(see lBrown. Hollander. & Korwarll974T:IWang & Wellsl2000l) 
with W r (2796) as the dependent variable. The null hypothesis 
of no correlation could not be ruled out for W,(2796) against 
L B /L* (0.8 cr), L K /L* K (0.2 a), M B (1.6 a), and M K (0.02 a). 
The lack of correlations found here are interesting in view of 
the compelling arguments by Bouche et al. (2006), who de- 
rive an anti-correlation between W, (2796) and galaxy lumi- 
nosity with a dependence on the faint-end slope of the galaxy 
luminosity function and a proportionality between galaxy lu- 
minosity and the cross-section of the absorbing CGM. 

To determine if W, (2796) has any dependency on galaxy 
color, we performed the non-parametric Kendall's r rank cor- 
relation test on W r (2796) against B-K. The test indicates that 
W r (2796) does not directly correlate with galaxy color for the 
full sample (Liu). A Kendall's r rank correlation test be- 
tween W r (2796) and B-K for the subsample with D < 50 kpc 
indicates that W r (2796) is not correlated with color at smaller 
impact parameters (0.3 a). 

The lack of a correlation be tween W r (2 7 96) an d B-K is 
consistent with the findings of IChen et al.l (l2010al) . who ex- 
amined B-R colors for 0.1 < z < 0.5. Howev er, the result 
is con trary to the statistically based results of Zibe tti etldl 
(120071) . who find bluer colors from the integrated light se- 
lected by stronger Mg II absorbers and redder colors selected 
by weaker absorbers for 0.4 < z < 1, where the minimum 
W, (2796) of their sample i s 0.8 A . A direct comparison with 
the work of IZibetti et all d2007l) is difficult due to its sta- 
tistical nature. However, a Kendall's r test examining the 
Wr(2796) > 0.8 A subsample of MAGllCAT yielded no cor- 
relation between W,(2796) and B-K (1.9 a). 

Bifurcating the full galaxy sample at W cut = 0.1, 0.3, 0.6, 
and 1.0 A, we find a 3.1 cr significance for a correlation be- 
tween W r (2796) and B-K for galaxies hosting W f (2796) > 
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1.0 A absorption; for the very strongest absorbers, the redder 
the galaxy, the greater the Mgll equiva lent width. As men - 
tioned, this is contrary to the findings of Zib etti et"aT1 ([2007). 

To examine the possible W f (2796) equivalent width dis- 
tribution dependence on B-K color, we also performed a 
Kolmogorov-Smirnov (KS) test comparing the cumulative 
W r (2796) distribution functions of red (B-K > 1 .39) and blue 
(B — K < 1.39) galaxies. We limited the test to only include 
those galaxies with detected Mg II absorption. The two sub- 
samples are statistically consistent with having been drawn 
from the same parent distribution (~ 0.5 cr). Limiting the 
sample to galaxies with D < 50 kpc, the KS test results remain 
below 1 (7. Though the distribution for red galaxies has power 
beyond the largest and smallest W r (2796) of the blue galax- 
ies, the statistics do not ferret out a difference between the 
two galaxy subsamples; the average and maximum W, (2796) 
for D < 50 kpc associated with red galaxies and blue galaxies 
are consistent (average 0.97 and 0.91 A, with maximum 2.9 
and 2.1 A, respectively). 

Our findings that the W r (2796) distributions from the CGM 
within D = 50 kpc of blue and red galaxies ar e indistinguish- 
able is at odds with the dramatic finding of iBordoloi et al.l 
(2011), who report a factor of eight larger equivalent width 
associated blue galaxies compared to red galaxies in stacked 
spectra for which the CGM is probed within D = 50 kpc. 

3.3. W,(2796) and Impact Parameter 

A commonly known property of Mg II galaxies is the anti- 
correlation between W r (2796) and impact parameter D, (e.g., 
Lanzetta & Bowed [19901: iBergeron & Bo ssel fT99lt ISteidell 
1995b IBouche et all 120061: IKacprzak et all l2008t IChen et al.l 
201 Oat IChurchill et al.l [2012a). In Figure Q] we present 
W r (2796) versus D. Galaxies with detected Mg II absorption 
are plotted as red points and those with upper limits on ab- 
sorption are plotted as blue points with downward arrows. 

We performed a non-parametric Kendall's r rank correla- 
tion test on W, (2796) against D, allowing for upper limits on 
W,(2796). For this sample, W,(2796) is anti-correlated with 
D at the 8.2 a level. Since W^(2796) cor relates with the num- 
ber of clouds (Voigt profile components. l Petitiean & Bergeron 
[19901 IChurchill et al.ll2003l lEvansll201 ll) . this result indicates 
that either the column densities, velocity spreads, or both, di- 
minish with projected distance from the galaxy. 

To parameterize the general behavior of the anti- 
correlation, we used the Expectat i on-M aximization 
maximum-likelihood method (Wolynetz 1979) account- 
ing for upper l imits on W r (2796) to fit a power law (see 
ICnenetal.lT2010al). logW r (2796) = ailog£>+a 2 , and a log- 
linear fit (see Churc hill et alj2012al) . logW f (2796) = aiD+a 2 , 
to the data. The data are poorly described by a power law 
due to the substantial population of W r (2796) < 0.1 A 
absorbers and upper limits. We present the log-linear fit and 
its 1 a uncertainties in Figure [T]for a\ = -0.019 ±0.002 and 
a 2 = 0.28 ±0.11. 

The considerable scatter about this relation may sug- 
gest that W,(2796) is governed by physica l processes re- 
lated t o the galaxy such as luminosity (cf., Kacprza k et al.l 
1 20081: IChen et all 1201 OaT). mass (cf.. IBouche et al.1 120061: 
iGauthier et al.l 120091: IChurchill et al.l 1201 2bl) . star formation 
(cf.. IChen et al.l l2010bt iMenard et al.l 1201 ll). or orientation 
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FIG. 1. — The Mgll A2796 rest-frame equivalent width, W,(219(>), ver- 
sus impact parameter, D. Galaxies with detected Mgll absorption are pre- 
sented as red points, whereas those with upper limits on absorption are blue 
points with downward arrows. The full sample comprises 169 galaxies with 
an 8.2 a anti-correlation between W,-(2796) and D. The solid line is a log- 
linear maximum likelihood fit to the data, logW r (2796) = ct[D + ct2, where 
Oi i =-0.019±0.002 and a 2 = 0.28±0.11. Dashed curves provide 1 cr un- 
certainties. 



(cf.. IBordoloi et al.ll201 It IBouche et al.[l201 lHKacprzak et al 
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Alternatively, the CGM is inherently patchy 



(IChurchill et alJl2Wl2^12ah . 

To examine possible trends and/or sources of this scatter 
in the distribution of W r (2796) versus D, we applied the two- 
dimensional Kolmogorov-Smirnov (2DKS) test. This test ex- 
amines if the 2D distribution of W, (2796) and D from two 
samples can be ruled out as being consistent. 

We first investigated bifurcations of the full sample about 
the median values of B-K, z ga i, L B /Lg, and L K /L* K , perform- 
ing a 2DKS test on the subsamples. Blue galaxies predom- 
inantly populate the region W r (2196) > 0.3 A and D < 40 
kpc, whereas red galaxies populate the W r -D plane more 
uniformly; however, the trend is not statistically significant 
(1.7 cr). Using a redshift cut, the data show a weak trend for 
larger W,-(2796) at fixed D for high z galaxies (2.3 cr). A 2DKS 
test on the luminosities shows that high Ls/Lg and Lk/L* k are 
found with larger W,-(2796) at fixed D than low luminosities 
at the 3.1 cr and 3.3 a level, respectively. 

We also divided the sample into quartiles by L B / L B , Lk / L* K , 
B-K, and z g - d \ and performed the 2DKS test on the lowest and 
highest quartiles. We obtained P(KS) = 0.012, 0.019, 0.36, 
and 0.36, for L B /L* B , L K /L* K , B-K, and z ga i, respectively. At 
best, there is a 2.5 a correlation that the distribution or scatter 
of W r (2796) versus D has a connection to L B jL* B , and a 2.35 a 
correlation it is connected to Lk/L* k . There is clearly no color 
or redshift dependence. 

3.4. Covering Fraction and Impact Parameter 

Here we examined the dependence of covering fraction di- 
rectly on impact parameter and galaxy luminosity, color, and 
redshift. We present directly observable quantities, making no 
assumptions with regard to the Mg II absorbing CGM density 
profile or geometry. 

3.4. 1 . Covering Fraction Within Fixed Impact Parameters 
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FIG. 2. — The covering fraction, /omax, inside D maJt for different W r (2796) bifurcations at W cu t. Shaded regions represent 1 <r uncertainties based upon binomial 
statistics. — (upper) The full sample of galaxies (solid line), the high luminosity galaxies (dashed line), and the low luminosity galaxies (dotted line), divided at 
the median B-band luminosity, Lb/L* b = 0.526. — (middle) The blue (B—K < 1.39) and red (B—K > 1.39) galaxy subsamples bifurcated at the median galaxy 
color. — (lower) The high (z > (z)) and low (z < (z)) redshift galaxy subsamples where (z) = 0.346. 



In order to examine the covering fraction within a fixed 
projected radial distance from the galaxies, we computed the 
quantity / 0mm = f(W > Wc Ut ,D < D max ), which we define as 
the fraction of absorbers with W-(2796) > W cu t inside a pro- 
jected separation of Z) max from the galaxy. 

In Figure|2] we plot f Dmm against D max in 10 kpc intervals. 
The uncertainties in / flmax are shown as shaded regions. These 
are the upper and lower limits on f Dm ?t , calculated u sing the 
formalism for binomial statistics (see|Gehrels 1986). Values 
of / D max at Omax < 10 kpc are not robust given the small number 
of galaxies within this impact parameter 5 . Due to the relative 
undersampling of galaxies at Z) max > 150 kpc and the cumu- 
lative nature of f Dm ^, the covering fraction will be virtually 
unchanging outside this impact parameter; thus the shape of 
/ Cm „ versus Z) max is less robust for Z) max > 150 kpc. 

In Figure [2] (upper panels), we present / 0mal for the full 
galaxy sample and for low and high luminosity galaxies. For 
each W cut subsample, the covering fraction of the Mgll ab- 
sorbing CGM inside D = Z) max decreases as D max is increased. 
At a given D max , there is a clear trend that as W cu t is increased, 
/omax decreases, indicating that the covering fraction within a 
fixed projected separation increases as the minimum absorp- 
tion threshold is lowered. The covering fraction exhibits no 

5 Of the three galaxies at D < 10 kpc, one is at low redshift where the an- 
gular separation from the quasar is much greater than the quasar seeing disk, 
and the other two were found at high er redshift following point spre ad func- 
tion subtraction of the quasar iSteidel, Dickinson, & Persson 1994), which 
was not performed in all surveys. 



luminosity dependence for W cut = 0.1, 0.3, and 0.6 A. How- 
ever, for W cut = 1 .0 A, f Dmiti is larger for high luminosity galax- 
ies at all D max than for low luminosity galaxies. This differ- 
ence is accentuated at Z) max < 50 kpc. 

In the middle panels of Figure |2] we present / 0m „ and its 
uncertainties versus Z) max for red and blue galaxies. The mean 
redshift for red galaxies is z = 0.36 over the range 0.1 < z < 
0.85, and for blue galaxies z = 0.44 over 0.1 <z< 1.02. For 
each Wcut subsample, the blue galaxies have higher covering 
fractions relative to red galaxies for D max > 50 kpc. Within 
50 kpc, the larger uncertainties make it difficult to distinguish 
any possible differences. 

In Figure [2] (lower panels), we present / 0mal and its uncer- 
tainties against Z) max for the high z and low z galaxy subsam- 
ples. For D max < 50 kpc, it is difficult to distinguish any 
redshift evolution in / Cmax due to the uncertainties. Beyond 
50 kpc, the data suggest that / Dm .„ may evolve such that at 
higher redshift, the covering fraction of Mg II absorbing gas 
is higher than at lower redshift. The differences in redshift 
may be connected to the fact that a larger proportion of galax- 
ie s at z > (z) have br ighter M B (see § 13.11 ). 

IChen et all d2010al) also examined the luminosity depen- 
dence of / Dnmx as a function of W cu t, where their D mdx is fixed 
at the luminosity scaled "gas radius", with j? gas cx: (Lg/Ll) 035 , 
assuming isotherma l density profile ([Tinker & Chen! [20081) 
and NFW profile ( Navar ro" et al.l 1 19961) models. As such, a 
direct comparison with our non-parameterized results is diffi- 
cult. Nonetheless, they found that /„„,„ has little to no depen- 
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FIG. 3. — The covering fraction profile, /< ), for the impact parameter bins < D < 25 kpc, 25 < D < 50 kpc, 50 < D < 100 kpc, and 100 < D < 200 kpc 
for different W, (2796) bifurcations at W cu t- The horizontal bars indicate the impact parameter bin width and the vertical bars are the 1 a binomial uncertainties. 
Data points are plotted at the mean impact parameter of the galaxies in the bin. — (upper) The full sample of galaxies (solid points), the high luminosity galaxies 
(open triangles), and the low luminosity galaxies (open squares). — (middle, lower) Subsamples and colors are the same as given in Figure[2] 



dence on galaxy B-band luminosity for W cut = 0.1, 0.3, and 
0.5 A, which is consistent with our result; for W cut = 1 .0 A, we 
found systematically larger / 0m „ for high luminosity galax- 
ies as compared to low luminosity galaxies, especially for 
£>max < 50 kpc. 

Comparing the fraction of absorbers to non-absorbers as a 
function of azimutha l angle between the quasar sight l i ne and 
galaxy major axis, iKacprzak. Churchill. & Nielsen! (2012) 
find / flm » ~ 70% for D max ~ 100 kpc and W cut = 0.1 A. Their 
result is consistent with our finding. 

3.4.2. Covering Fraction Profiles 

To examine the covering fraction profile with projected dis- 
tance from the galaxy, we computed f {D) = f(W > W cut , (D)), 
which we define as the fraction of absorbers with W,(2796) > 
Wcut in fixed impact parameter bins. We select the bins 
< D < 25 kpc, 25 < D < 50 kpc, 50 < D < 100 kpc, and 
100 < D < 200 kpc, for which (D) is the average impact pa- 
rameter of the galaxies in the bin. 

In Figure [3] we plot f {D) against (D). The horizontal bars 
are the impact parameter ranges and the data points are the 
average impact parameters, (D), for each bin. The vertical 
error bars are the 1 a uncertainties in f {D) based upon bino- 
mial statistics. The subsamples presented in each panel are 
identical to the corresponding panels of Figure|2] 

As can be seen in the upper panels of Figure [3] the f (D ) 
profile decreases as impact parameter is increased and this 
behavior is exhibited regardless of W cu t- Also note that as 



TABLE 2 

Luminosity Dependence of Covering Fraction Profiles, 
Wcut, A (0-25), kpc [25-50), kpc [50-100), kpc [100-200), kpc 
All Galaxies 
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W cut is increased, f {D) is smaller for a given impact parame- 
ter bin. This indicates that in a fixed annulus around galaxies, 
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the sky-projected distribution of the Mg II CGM becomes pro- 
gressively patchier as higher column density material and/or 
more complex kinematics are selected. For W cu t = 0.1 A, the 
f {D) profile decreases from unity within a radius of 25 kpc to 
20% in the annulus 100-200 kpc, whereas for W cut = 1.0 A, 
the covering fraction decreases from 40% to less than 10%. 

Also shown in the upper panels of Figure [3] is f (D) split 
by Lb /L^ = 0.526. For low luminosity galaxies, the covering 
fraction vanishes outside of D = 100 kpc for all Wcut- On the 
other hand, in the annulus 100-200 kpc, the covering fractions 
for high luminosity galaxies are approximately 40%, 30%, 
10%, and 6% for W cut = 0.1, 0.3, 0.6, and 1.0 A, respectively. 
Covering fraction profiles for all, low, and, high luminosity 
galaxies are listed in Table [2] 

The data are suggestive of a luminosity dependence for 
the maximum extent of Mgll absorbing gas, such that the 
CGM of low Lb/L% galaxies extends no further than 100 kpc, 
whereas high luminosity galaxies have a detectable Mgll 
CGM beyond 100 kpc. The difference between high and low 
luminosity galaxies is most prominent outside of 50 kpc for 
Wcut = 0.1, and 0.3 A but is most prominent inside 50 kpc for 
Wcut = 0.6, and 1.0 A. As such, high luminosity galaxies ex- 
hibit a more extended Mg II CGM, whereas for low luminosity 
galaxies, the Mg II CGM is more concentrated in the central 
regions but with relatively lower covering fraction as W cu t is 
increased. 

In the middle panels of Figure [3] the covering fraction pro- 
files of red and blue galaxies exhibit no significant differ- 
ences. However, note that for red galaxies, none have ab- 
sorption with W r (2796) > 0.6 A in the range 100-200 kpc. 
Blue galaxies consistently have covering fractions of approx- 
imately 20% for all W C ut in the range 100-200 kpc. How- 
ever, within 25 kpc the covering fractions for blue galaxies 
drops more rapidly than for red galaxies with increasing W cut . 
Moreover, for W cut = 0.6 and 1.0 A there is a weak trend for 
the covering fraction of red galaxies to drop more rapidly for 
D > 25 kpc relative to the D = 0-25 kpc region. This may sug- 
gest that optically thick gas is relatively less abundant in the 
outer regions of redder galaxies as compared to the relative 
abundance in the outer regions of bluer galaxies. 

As shown in the lower panels of Figure [3] high redshift 
galaxies exhibit a trend of higher / (0> in the 25-50, 50- 
100, and 100-200 kpc annuli relative to low redshift galax- 
ies. Within 25 kpc, the covering fraction may be higher in 
low redshift galaxies for stronger absorption, W cut = 0.6 and 
1.0 A. Though the trends are marginal, they may suggest a 
greater extension of the gas at higher redshift followed by a 
settling of material into the inner regions at lower redshift. 
However, the behavior of f {D) with redshift may be partially 
influenced by the larger fraction of higher luminosity galaxies 
at z>(z). 

IChen et ail (|2010a) examined f (D) for their sample ((z) = 
0.25), and also found that the covering fraction profile, f (D) , 
decreases with increasing D and is smaller at a given (D) 
with increasing W cllt using W cllt = 0.1, 0.3, and 0.5 A. How- 
ever, a direct co mparison with our results is difficult because 
IChen et ail (1201 Oal) scale D assuming a Z?-band galaxy lumi- 
nosity dependent "gas radius" proportional to (Lb /Lg) 035 . 

3.5. Luminosity Scaling and Covering Fraction 

Our results in the previous section strongly suggest that the 
Mg II absorbing CGM extends further for weaker absorption 
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Luminosity Scaled Halo Absorption Radii 
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than for stronger absorption, and that this behavior has a clear 
dependence with galaxy luminosity. What is further clear, 
is that for a given W C ut, the CGM becomes more "patchy" 
with increasing D, as indicated by the decreasing covering 
fractions. As such, it would seem that the notion of a well- 
defined Mgll CGM "gas radius", or "halo absorption ra- 
dius" is an oversimplification in that the extent of the gas ex- 
hibits a "fuzzy" boundary when averaged over many galax- 
ies. Nonetheless, we parameterize the characteristics of an 
"outer" boundary and its plausible dependence on the galaxy 
properties and redshift 

Since the work of .Bergeron & Boisse (1199 lh . the extent of 
absorbing gas is commonly assumed to follow a Holmberg- 
like relation, R(L) = R«(L/L*)P. We examined whether the 
halo absorption radius also depends on W-(2796), and/or 
galaxy color and redshift. 

To examine the dependence on W r (2796), we adopt W cu t = 
0.1, 0.3, 0.6, and 1.0 A. To estimate and (3, we varied the 
two parameters over the ranges < /?* < 300 kpc and < j3 < 
1 and computed the function q(R*,(3,L, W cut ) = wr(>) + r(<), 
where r ( >) is the fraction of systems with W,(2796) > W cu t 
below R(L), r(<) is the fraction with W-(2796) < W cu t above 
R(L), and w is a weighting factor 6 . The parameter values 
are adopted when q(R*,(3,L, W cut ) is a maximum, for which 
the covering fraction inside R(L) is / S(i) = 1/(1 +;t), where 
x= («(>)/«(<))(l -r(<))/V(>), and where «(>) and «(<) are the 
number of systems with W f (2796) > W cu 7 and W r (2796) < 
Wcut, respectively. 

The downward and upward 1 a uncertainties in /?» were 
estimated by using the best estimate of f3 and performing 
one-sided integration under the q(R*,(3,L, W cu t) curve until 
84.13% of the total area was obtained. Similarly, the down- 
ward and upward 1 a uncertainties in (3 were estimated by 
using the best estimate of The 1 a uncertainties in f m) are 
computed using binomial statistics by computing the fraction 
of galaxies with W, (2796) > W cut to all galaxies within R(L). 

Our results are presented in Table[3]for both B- and K-bwA 
luminosities. In the upper panels of Figure |H we present D 
against Lb)L* b . The luminosity scaling remains roughly con- 
stant at (3 ~ 0.4 for increasing W cu t. The absorbing gas halo 
radius, for an L* B galaxy is on the order of ~ 75 kpc for 
all but Wcut = 0.1 A, which is on the order of ~ 90 kpc. The 
covering fraction decreases from f KL) = 0.89 for W cu t = 0.1 A 
to f m = 0.34 for Wcut = 1 .0 A. 

In the lower panels of Figure 21 we present the Zf-band re- 
sults. We note that this subsample has 20 fewer galaxies than 
the fi-band subsample. The luminosity scaling flattens from 

6 Since we wish to determine the "outer envelope" of the halo absorption 
radius, we adopt the weight w -2. When w = 1, we find (3 consistent with 
zero. 




FIG. 4. — Impact parameter, D, versus luminosity for different W,-(2796) bifurcations at W cut . Green points are galaxies with W r (2796) < W cllt and purple 
points are galaxies with W r (2796) > W cut . — (upper) B-band luminosity. — (lower) X'-band luminosity. Open points represent galaxies which have limits on 
W,-(2796). The fit parameters R t , ($, and f RiL) give the absorbing gas halo radius for an L* galaxy, luminosity scaling power index, and absorption covering 
fraction, respectively, for each W r (2796) = W cu t bifurcation (see Table|5J. Dashed lines provide the envelope 1 cr uncertainties in the fit at fixed L/L* . 



(3 ~ 0.25 for W cut = 0.1 and 0.3 A to (3 ~ 0.15 for W cut = 0.6 
and 1.0 A, whereas B* decreases from ~ 75 kpc to ~ 60 kpc. 
The covering fraction behaves similarly to the B-band (see 
T ableEDl 

ISteidell {[995) reported /3 ~ 0.2 and B» = 55 kpc for the B- 
band and (3 ~ 0.15 and B» = 58 kpc for the /f-band, both for 
W cut = 0.3 A (the values of B* quoted here have been converted 
to a "737" ACDM cosmology at the mean redshift, (z) = 0.65, 
of his sample). He deduced a covering fraction of f m ~ 1 
and that (3 an d R* for the AT-band are slightl y smaller than for 
the B-band. iGuillemin & Bergeronl (11997b obtained similar 
results, with B* = 47 kpc (adjusted from qo = 0.05 and h^o at 
z = 0.7) and (3 = 0.28 for the B-band with W C ut ^ 0.3 A. 

Using the observed Mg H absorber red s hift p ath num- 
ber density as a constraint, Kacp rzak et al.l d2008l) explored 
the parameter space of the minimum luminosity of Mgll 
absorption-selected galaxies, the slope, (3, the halo absorption 
radius, and covering fraction, f m . For W cu t = 0.3 A, they 
constrained R* to the range 50-1 10 kpc, (3 to the range 0.2- 
0.28, and f R(L) to the range 50-80%. Our values of B* reside 
in the range of allowed values from their study. 



IChen et all (1201 Oal) reported fi* = 107 kpc with f3 = 0.35 
for the B-band, where they fitte d their data in the con text of 
an isothermal sphere model (see lTinker & Ch en 2008) of the 
Mgll CGM. For these values, they find f R{L) = 0.7 for W cu t = 
0.3 A and / fl(il = 0.8 for Wcut = 0.1 A, both of which are - 10% 
smaller than what we find. 

Emp loying the | Tinker & Chen! d2008l) isothermal sphere 
model, Bordol oTet al.l (1201 ll) deduced R* = 115 kpc. Whereas 
(3 for this work i s consistent with the value reported by 



Chen et al 



Chen et al 



(2010a), the larger values of R* re ported by both 
(1201 Oal) and lBordoloi et al.ld201 lb may be an ar- 
tifact of their applic ation of a model to th e data. Further- 
more, the methods of iBordoloi et al.l (1201 ll) involve "averag- 
ing" over annuli of fixed impact parameter, which may intro- 
duce an additional systematic. 

3.5.1. Galaxy Color 

To examine the dependence of R(L) on galaxy color, we 
computed B», (3, and f R(L) for red and blue galaxies, bifur- 
cated by B-K = 1 .39, in both the B- and ZT-bands. Due to the 
smaller number of galaxies in the blue and red subsamples (74 
each), we adopted a single equivalent width cut, W cu t = 0.3 A 
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FIG. 5. — Impact parameter, D, versus Lg/Lg (top) and LnjL* K (bottom) 
for blue and red galaxies bifurcated by the median color B—K = 1.39. Point 
types and curves are the same as Figure gjfor W cul = 0.3 A. The luminosity 
scaling, f3, and radius, R„, of the absorbing gas halo radius exhibits a color 
dependence for this W cll t; blue galaxies have a steeper luminosity scaling, 
larger R t , and larger covering fraction than red galaxies. 



(the median W,-(2796) for galaxies with measured colors). As 
shown in Figure[5]and tabulated in Table|4j the parameters for 
R(L) suggest a dependence on galaxy color in both the B- and 
/T-bands. 

In the B-band (top panels of Figure [5]), we find that blue 
galaxies have a steeper luminosity scaling, (3 ~ 0.45, and 
larger covering fraction, f m ~ 0.9, than red galaxies, (3 ~ 
0.3 and f mL) ~ 0.75. The halo absorption radius for blue 
galaxies, ~ 90 kpc, is larger than for red galaxies, ~ 
75 kpc. However, both values of R* are consistent within 
uncertainites. These values are smaller than those found by 
iBordoloi et alj d201 lb . who report 107 kpc for blue galaxies 
and 118 kpc for red galaxies. Furthermore, we find the oppo- 
site behavior of R* with color in that our blue galaxies have 
larger R* than red galaxies. 

The R{L) dependence on galaxy color is similar in the K- 
band (bottom panels of Figure |5]). The halo absorption ra- 
dius, i?*, is larger for red galaxies (though the values are 
still consistent within uncertainties), with R* ~ 80 kpc and 
~ 55 kpc for blue and red galaxies, respectively. The lu- 



FlG. 6. — Impact parameter, D, versus Lg/Lg (top) andLjc/L^ (bottom) for 
low and high redshift galaxies bifurcated by the median redshift z gll { = 0.346. 
Point types and curves are the same as Figuref4]for W cut = 0.3 A. The average 
Z for the low redshift subsample is 0.22 and the average for the high redshift 
subsample is 0.61, which is a ~ 3.2 Gyr difference. 

minosity scaling is steeper in blue galaxies, f3 ~ 0.3, than red 
galaxies, j3 ~ 0.1, and the covering fraction is larger in blue 
galaxies than red galaxies (f m ~ 0-9 and / fi(i) ~ 0.8, respec- 
tively). 

3.5.2. Galaxy Redshift 

We also examined the R(L) dependence in both the B- and 
/iT-bands on galaxy redshift, bifurcating the sample into low 
and high redshift galaxies by the median redshift z = 0.346, 
and adopting W cut = 0.3 A. Values for R*, (3, and f m are pre- 
sented in Figure|6]and tabulated in Table|4] 

In the B-band (top panels of Figure|6]l, we find R« ~ 90 kpc, 
(3 ~ 0.5, and f m ~ 0.8 in the low z subsample. The high 
z subsample has similar values for R* (~ 85 kpc) and f m 
(~ 0.85), but a much shallower dependence on Lb/L^, with 
f3 ~ 0.3. This difference in j3 is driven by the presence of Lg ~ 
Q.\-QAL* B galaxies with D ~ 30-50 kpc at higher redshift, 
which are less abundant in the low redshift subsample. It is 
doubtful that the possible evolution in (3 is a selection effect, 
since lower luminosity galaxies would presumably be selected 
against at high redshift, not lower redshift. 
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TABLE 4 

Luminosity Scaling for Subsamples at Wait = 0.3 A 
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- /f-band 




fan 


All 
Blue 
Red 

Low z 
High z 


79+50 

75 +57 
90 + J 5 6 

85+f 1 


n -is+o.oi 

0.45^ 
0.28^ 

53 +0 01 

n 97+O.O6 
U,z '-O.03 


0.84^ 

7S+0.0V 
"•' J -0.08 

n a i+o os 

U - S1 -O.07 

a os+0 05 


7s;f 

83!f 

55+ 54 
-"-17 

78!f 


0.28^$ 

0-H^oi 
0.2JC2S 


0.83^ 

r. Q9+0.04 
u -* z -0.06 

0.80^ 
0-83^ 



On the other hand, the /T-band (bottom panels of Figure [6]l 
shows no discernable difference between the low and high z 
subsamples. The values for (3, and / s(i) for both subsam- 
ples are ~ 80 kpc, ~ 0.3, and ~ 0.8, respectively. Compar- 
ing to the apparent Z?-band evolution, this may imply that the 
luminosity dependence of the extent of the Mgll absorbing 
CGM is unchanged with redshift for the K-b&nd. 

3.6. Luminosity Functions 

Prior to measurements of the l uminosity functions, 4>(M), 
of ga l axies out to z = 1 (e.g.. iLillv et al.l 119951: iLin et al.l 
19991 iFried et"H]l200Tt IWolf et al.l 12003b iFaber et al.l 12007b 
Cirasuolo et al.l 1201 Oh . selecting galaxies by Mg II absorp- 
tion provided a compelling technique for compiling a pre- 
sumably complete sample of interm ediate redshift galaxies 
(ISteidel. DickinsonT& Persson 1994L Using absorption se- 
lection, ISteidel. Dickinson. & Perssonl (1 19941) presented the 
B- and K-b&nd luminosity functions of intermediate redshift 
galaxies associated with W,(2796) > 0.3 A. 

In order to compare to and ex pand upon the work of 
ISteidel. Dickinson. & Pers son (1994), we determined the B- 
and A'-band luminosity functions for MAGllCAT. We divided 
the galaxies into six subsamples, bifurcated by W r (2196) = 
0.3 A for all redshifts, low z, and high z- The average red- 
shift for the low z subsample is z = 0.22, and for high z is 
z = 0.61, which translates to a 3.2 Gy r time spread. Following 
ISteidel. Dickinson. & Perssonl d 19941) . we applied a 1 /R 2 (L) 
correction factor to account for the dependence of gas halo 
radius on galaxy luminosity, using f3 = 0.38 for the fi-band, 
and j3 = 0.28 for the /T-band, both from Table[3] This rectifies 
the relative volume probed by absorption line surveys at fixed 
luminosity under the assumption of completeness. 

In Figures|7]and|8] we present the B- and /^-band luminosity 
functions. For each binned data point, we have computed the 
mean B-K color of the galaxies contributing to the bin. Blue 
data points indicate a mean blue color (B-K < 1 .39) and red 
points indicate that the mean color is red (B-K > 1.39). For 
reference, we have overplotted the Schechter luminosity func- 
tions (plus a single additive constant to roughly match the data 
at M* ) for both z = 0-3 (s olid curve) and z = 1 . 1 (dashed curve) 
from Faber et ail d2007l) for the Z?-band and ICirasuolo et al.l 
(2010) for the Zf-band. If the galaxy surveys are magnitude- 
limited, we estimate that the completeness of MAGllCAT be- 
gins to decline for Mb > -18 and Mk > -17.5 for the low z 
subsample and Mb > -19 and Mk > -20 for the high z sub- 
sample; these data are conservatively plotted as open points. 

3.6.1. Comparing to Steidel et al. (1994); W, (2796) > 0.3 A 

The W r (2796) > 0.3 A subsample at high z (Fig- 
ure |7]f) is most appropriately compared to the 



ISteidel. Dickinson. & Perssonl (1 19941) sample. General 
trends noticeable for galaxies in the B-band associated with 
W r (2796) > 0.3 A absorption are that there is a slightly 
higher relative abundance of luminous galaxies in the high 
z subsample compared to low z (Figures Wf-d), whereas 
in the /T-band, the most luminous galaxies (Mk < -24.5) 
associated with W r (2796) > 0.3 A are not selected in both 
low and high z subsamples (Figures |8}/-/). This is consistent 
with the correlation between Mb and redshift we found (see 
§ 13. U . For all subsamples, it also appears that the faint-end 
slope in the B-band may be slig htly flatter than the luminosity 
functions of IFaber et al.l (120071) . In fact, the measured <&(Mb) 
for W f (2796) > 0.3 A at high z (Figure |7]f ) is suggestive 
of the turn down at Mb > —20, w hich was reported by 
ISteidel. Dickinson. & Perssonl (11994 . 

Apart from the lack of Mk < -24.5 galaxies for the 
W r (2796) > 0.3 A subsamples, the /T-band luminosity func- 
tions a ppear to more or less match the luminosity func- 
tions of Cirasuolo et al. (2010). We see a possible turnover 
in the faint-end slope for the W r (2196) > 0.3 A, high 
z subsample (Figure [HQ. This result is contrary to 
ISteidel. Dickinson, & Perssonl (1 19941) . 

For the fi-band luminosities, we note that the observed 
trends (i.e., flattening of the faint-end slope and rela- 
tive overabundance at high luminosity) are reminiscent of 
the Malmquist bias that plagues magnitude-limited surveys. 
Given the heterogeneous selection methods used by the vari- 
ous works from which we constructed MAGllCAT, it is diffi- 
cult to quantify the degree to which this may be an issue. 

3.6.2. Redshift and W r (2796) Differences 

We ran KS tests on the unbinned data to investigate pos- 
sible statistical differences between the observed luminos- 
ity functions for each subsample. For the Z?-band, both the 
W r (2796) < 0.3 A and W r (2796) > 0.3 A subsamples are 
consistent with being drawn from the same population, with 
P(KS) = 0.27 for all redshifts (Figures|7b-fc), P(KS) = 0.19 for 
low z (Figures|7]:-(i), and P(KS) = 0.42 for high z (Figures|7^- 
/). The same tests in the /^-band yielded similar results for 
the full sample with P(KS) = 0.41 for all redshifts (Figures|HJj- 
b), P(KS) = 0.66 for low z (Figures ®-d), and P(KS) = 0.74 
for high z (Figures^-/). 

Comparing the low z and high z subsamples (with mean 
redshifts of 0.22 and 0.61, respectively, which is a ~ 3.2 
Gyr difference), we found evidence for redshift evolution 
only in the B-band, $(Mb), for the galaxies associated with 
W,(2796) < 0.3 A absorption (Figures |7J;-e); the test yielded 
P(KS) = 0.0003 (3.6 a). This result is due to the absence of 
luminous galaxies with Mb < -22.5 and the relative paucity 
of galaxies with -22.5 < Mb < -21 at low z, and a lack of 
lower luminosity galaxies with Mb > -18 at high z- The sig- 
nificance level of the KS test may be accentuated by possible 
magnitude limiting effects in the high z subsample (but the 
bright end would not be related to this plausible selection ef- 
fect). 

The same trend is also apparent for the W,(2796) > 0.3 A 
subsample in the B-band (Figures |7J/-/). However, the re- 
sult is not formally significant; the KS test yielded a 91% CL 
(1.7 a). For $(M^), there is no suggestion of redshift evolu- 
tion for both the W,.(2796) < 0.3 and W r (2796) > 0.3 subsam- 
ples. 

3.6.3. Color Sequence Along <3>(M B ) and $(M K ) 
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FIG. 7. — B-band luminosity functions for W r (2796) < 0.3 (a, c, and e) and 
W r (2796) > 0.3 (i, rf, and /). Galaxies at all redshifts are included in panels 
a and b, while galaxies with z < (z) are in panels c and d (mean z = 0.22), and 
galaxies with z. > z) are in panels e and / (mean z = 0.61), where (z) = 0.346. 
The solid and dashed curves are z = 0.3, 1.1 Schechter functions, respectively, 
from Faber et al. 1 2007). Data points are colored according to the mean color 
of the galaxies in the bin. Blue points correspond to a mean color B—K< 1 .39 
and red data points correspond to B — K > 1.39. Open points are where the 
completeness of the sample declines. 



For the B-band, it is somewhat striking that, for high z 
galaxies associated with W r (2796) < 0.3 A (Figure |7^), the 
mean color trends toward bluer galaxies for a wide range 
of luminosities. Yet, at low z, the mean color of the galax- 
ies associated with W r (2796) < 0.3 A absorption (Figure [TJr) 
trends toward redder galaxies. For W, (2796) > 0.3 A, galaxies 
fainter than M B trend red at lower z (Figure|7J/). However, the 
color sequence along $(M B ) for the low redshift and high red- 
shift subsamples (Figures |7J/ and [7)0 show no clear redshift 
tr end. 

ISteidel. D ickinson. & Perssor] d 19941) reported a correlation 
between Mk and B-K such that lower luminosity galaxies in 
the K-b&nd tend to be more blue. Performing a Kendall's r 
rank correlation test on Mk versus B-K for the full sample 
of galaxies presented here yielded a 6.3 a correlation. For the 
^-band, the trend for lower luminosity galaxies to be bluer 
is clearly seen in the color sequence of $(Mk) for all red- 
shift and W r (2796) subsamples (consistent with the correla- 
tion tests on Mk versus B-K). There is no correlation be- 
tween B - K and M B . 

To further explore the data, we bifurcated the sample at 
(1) W cut = 0.3 A, (2) z = 0.346, and (3) L B /L* B = 0.526. The 
Kendall tests on M K versus B-K for these subsamples re- 
sulted in significant correlations in all cases, with the most 
significant results for low L B jL* B (5.9 a), high z (5.3 a), and 
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FIG. 8. — ^-band luminosity functions for W r (2796) < 0.3 (a, c, and e) 
and Wr(2796) > 0.3 (b, d, and /). Galaxies at all redshifts are included 
in panels a and b, while galaxies with z < (z) are in panels c and d (mean 
z = 0.22), and galaxies with z > (z) are in panels e and / (mean z = 0.61), 
where (z) = 0.346. The sol id and dashed curves are z = 0.3, 1.1 Schechter 
functions, respectively, from Cirasuolo et al. ( 2010). Point types are the same 
as described in Figure[7j 



W,-(2796) > 0.3 (4.6 a). These subsamples all correspond to 
the faint end of &(Mk) shown in Figure[8]f . 

4. DISCUSSION 

With the larger sample afforded by MAGllCAT, we have 
obtained additional leverage to probe the dependence of the 
properties of the Mg II absorbing CGM with galaxy color, lu- 
minosity, impact parameter, and W r (2796) threshold. For our 
analysis, we have purposely avoided couching the data within 
the framework of model expectations or scaling various mea- 
sured quantities with "second parameters", such as scaling the 
impact parameter by galaxy luminosity, in order to provide a 
direct view of the CGM-galaxy connection and to interpret 
the data from a purely observational perspective. 

With regard to the relationship between galaxy color and 
W,-(2796), our results challenge those of previous studies. For 
the overall sample, we find no correlation between W r (2796) 
and B-K. However, we find that for the subsample with 
WX2796) > 1 A, redder galaxies have larger W r (2796) with 
a significanc e level of 3.1 a. Th ese findings are contrary t o 
the results of IZibetti etail d2007l) and lBordoloi etalj d2011l) , 
who both report larger W, (2796) associated with bluer galax- 
ies. It is difficult to rectify the inconsistencies of our findings 
with those works unless the solution originates in the different 
e xperimental method s . 

Zibe tti et"aT] ((2007) developed an image stacking method in 
which the individual galaxies are not identified, but light from 
galaxies in various SDSS photometric bands are measured rel- 
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ative to control fields for which the quasars show no Mg II 
absorption. They employ observer-frame colors, whereas we 
examine rest-frame colors in cluding the K-bsnA. The color 
cuts of the galaxies studied by Bord oloi et al.l(l201 ll) are based 
upon a color-mass relation using u -B colors. Furtherm ore, 
the equivalent widths measured by iBordoloi et al.l (1201 ll) are 
based on stacked spectra over fixed impact parameter annuli 
in which the Mg II doublet is not resolved, whereas those mea- 
sured in this work are based upon pencil-beam probes of the 
CGM for which all double ts are cleanly resolved. If the CGM 
is intrinsically patchy (see Churchill et al. 2012a) with cover- 
ing fraction dependencies on W r (2796) threshold and galaxy 
luminosity (as we have demonstrated in this paper), then the 
two methods (statistical versus case-by-case) may be provid- 
ing complimentary but different clues to the nature of the 
CGM. 

Our analysis of the covering fraction profile of the Mgll 
absorbing CGM as a function of impact parameter, galaxy 
luminosity, and absorption strength threshold provides im- 
proved insights on the extent and distribution of the low ion- 
ization, metal-enriched CGM (see Figure[3]and Table|2]i. The 
decrease in the covering fraction with increasing absorption 
threshold indicates that the cross-section of the highest col- 
umn density and/or most kinematically dispersed gas (plau- 
sibly winds material) is smaller than the cumulative cross- 
section including lower column density, kinematically quies- 
cent gas (plausibly accretion material). 

Additional insight is provided by the relative behavior of 
the covering fraction profile of low and high luminosity galax- 
ies (bifurcated at L B /L* B = 0.526). First, the Mgll absorbing 
CGM appears to extend no further than 100 kpc for low lumi- 
nosity galaxies. Second, the CGM at D > 100 kpc surround- 
ing high luminosity galaxies is dominated by lower column 
density, kinematically quiescent gas traced by W,(2796) < 
0.6 A absorption; the cross-section of this material increases 
as the equivalent width threshold is decreased to 0. 1 A. This 
trend is also seen in the range 50 < D < 100 kpc. In con- 
trast, in the range D < 50 kpc, the observed frequencies of 
lower and higher column density CGM gas is not strongly 
dependent on galaxy luminosity for W r (2796) < 0.6 A, but 
for W r (2796) > 0.6 A, high luminosity galaxies have a higher 
observed frequency of Mgll absorbing CGM. This differ- 
ence is most highly apparent for W r (2796) > 1.0 A. Within 
D = 25 kpc, the covering fraction is effectively unity for 
the thresholds Wr(2796) o > 0.1 and 0.3 A, but declines for 
W, (2796) > 0.6 and 1.0 A as the threshold is increased. 

These results indicate a strong luminosity dependence in 
cross-sections and extent of the Mgll absorbing CGM such 
that high luminosity galaxies have both a much more extended 
diffuse CGM than low luminosity galaxies and higher filling 
factors of the highest column density material within 50 kpc 
as compared to low luminosity galaxies. Since the luminos- 
ity segregation we are employing is for the B-band, we could 
be seeing a connection between the CGM and young stars, 
such that the increased presence of young stars is related to a 
greater cross-section of higher column density gas in the in- 
ner regions of the galaxies and a more extended distribution 
of the lower column density gas in the outer regions of the 
galaxies (beyond 100 kpc). This would imply a stellar driven 
mechanism for distributing the gas in the CGM. 

Following previous studies, we have assumed the relation 
R(L) = R*(L/L*Y to describe the "halo absorption radius" de- 
pendence on luminosity. For the B-band luminosity, we find 



that the sensitivity of R(L) to luminosity (the parameter f3, see 
Tableland Figure @]i does not change as the equivalent width 
threshold is raised. For W r i2196) > 0.1 we find /3 ~ 0.36, 
and for W r {2191) > 0.3, 0.6, and 1.0 A, we find (3 ~ 0.38. 
However, the extent for an L* B galaxy tends to decrease as the 
equivalent width threshold is raised. We find B* ~ 90 kpc for 
W r (2796) >0.1A,S»~ 80 kpc for W r (2796) > 0.3 A, and 
B» ~ 75 kpc for W r (2196) > 0.6 and 1.0 A. 

Interestingly, we find a shallower dependence of R(L) on 
the /f-band luminosity. Contrary to the B-band, the value of 
(3 tends to decrease as the equivalent width threshold is in- 
creased (though not monotonically). Furthermore, we find 
smaller B* for the /f-band than for the B-band. 

The greater extent and higher luminosity sensitivity to the 
B-band relative to the Zf-band would further strengthen the 
idea that the geometric extent and morphology of the Mg II ab- 
sorbing CGM is governed in large part by young stars. There 
is an indication that galaxy color directly plays a role in gov- 
erning the higher B-band luminosity sensitivity of R{L). For 
W r {2191) > 0.3 A, blue galaxies have a remarkably steeper lu- 
minosity dependence than do red galaxies. For blue galaxies, 
we find (3 ~ 0.45 for the B-band and f3 ~ 0.30 for the /T-band, 
whereas for red galaxies, we find f3 ~ 0.28 for the B-band and 
(3 ~ 0.11 for the Zf-band. In the B-band, the covering frac- 
tion within R(L) of red galaxies is ~ 15% lower than for blue 
galaxies, whereas it is only ~ 10% lower in the Zf-band. 

We strongly caution that the function R(L) should not be 
viewed as a well-defined outer boundary to the Mgll ab- 
sorbing CGM. Examination of the data on the D-L/L* plane 
shown in Figure|4]clearly reveals that a non-negligible number 
of galaxies with W, (2796) greater than the threshold cut also 
reside above the curve. The 1 a upper envelope of the R(L) 
boundary characterizes the presence of these galaxies. These 
galaxies also drive the large upward uncertainties in best fit 
B*. Note that this uncertainty increases from ~ 36 to ~ 57 kpc 
as the absorption threshold is decreased from W r (2796) > 1 .0 
to 0.1 A. This suggests that the R(L) boundary is less well 
defined for the lower column density structures of the CGM. 
Physically, the behavior of the data in this regime of impact 
parameter and W r (2196) could either reflect a greater level of 
patchiness in the low column density material residing in the 
outskirts of the CGM of individual galaxies or a broad range 
of CGM properties from galaxy to galaxy. 

The color sequence across the B- and Zf-band luminosity 
functions indicates relationships between W,-(2796), redshift, 
and galaxy luminosity that might otherwise not have been 
discovered in the data without a multivariate analysis. At 
high z, W,(2796) < 0.3 A absorption is associated primarily 
with galaxies for which the mean color is blue (B-K < 1 .39) 
across the full B-band luminosity range. Interestingly, we 
find the mean colors of the low z galaxies associated with 
W,-(2796) < 0.3 A absorption have an average color that is 
primarily red. For W r (2796) > 0.3 A, the data may suggest 
that the reddening of galaxies with decreasing redshift occurs 
at sub-Lg luminosities. 

Comparing to the luminosity functions of 
ISteidel. Dic kinson. & Perssonl (Il994h . we see possible 
evidence for a faint-slope turnover in the B-band for the 
appropriate subsample of galaxies in MAGllCAT (see Fig- 
ure UY)- Interestingly, we also see possible evidence for a 
faint-slope turnover in the corresponding Zf-band luminosity 
function at the faint end (see Figure [8]/'). This latter result 
is contrary to the findings of Steidel, Dickinson. & Perssonl 
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(1994). Given the heterogeneous selection of the galaxies 
for MAGllCAT, it is difficult to assess if these turnovers are 
influenced by incompleteness at low luminosities at higher 
redshift. 

5. CONCLUSION 

Combining our previous studies and the extant works in the 
literature, we have compiled a "Mg II Absorber-Galaxy Cata- 
log" (MAGllCAT) of intermediate redshift galaxies and their 
associated circumgalactic medium (CGM) as probed using 
Mgll AA2796,2803 absorption. Details of the MAGllCAT 
data will be presented in Paper I. In this paper, we present re- 
sults from a first analysis, focusing exclusively on "isolated" 
galaxies, which are defined to have no neighboring galaxy 
within a projected distance of 100 kpc and line of sight ve- 
locity within 500 km s" 1 . 

The sample presented here comprises 169 galaxies toward 
117 quasar sightlines over the redshift range 0.072 < z ga i < 
1.120, with median (z) = 0.346. The rest-frame magnitudes 
range from -16.1 > M B > -22.8 (0.017 < L B /L* < 7.44) and 
-17.0 > M K > -24.3 (0.006 < L K /L* K < 4.9) with AB rest- 
frame colors 0.04 < B-K < 4.09. The median B-band lumi- 
nosity is Lg /L b = 0.526 and the median color is B-K = 1.39. 

The main results are: 

1. The B-K color exhibits a trend to redden with decreas- 
ing redshift, but the reddening is not statistically significant 
(2.3 ct). The mean Mb increases with increasing redshift 
(4.4 a), whereas Mk exhibits a weak trend to increase with 
redshift (2.1 a). The subsample with W r (2796) > 0.8 A is 
highly c onsistent with no col or evolution, in conflict with the 
result of Zibett i et al.l d2007l) : however, the subsample with 
W,-(2796) < 0.8 A shows a reddening trend with decreasing 
redshift (2.5 a), suggesting the reddening trend we see is 
dominated by galax ies with W r (2796) below the sensitivity of 
IZibettietalJ (120071) . 

2. There is no correlation between W r (2796) and B-K 
for the full sample (1.1 a). However, for galaxies asso- 
ciated with W,-(2796) > 1.0 A absorption, B-K correlates 
with W r (2796) at 3.1 a, indicating that redder galaxies have 
stronger absorption in this equivalent width regime. In our 
sample, the distributions of W, (2796) within D = 50 kpc of 
blue (B-K < 1.39) and red (B-K > 1.39) galaxies are con- 
sistent with being drawn from the same parent distribution 
(0.5 a). The W r (2796)-color correlation in our sample con- 
flicts with the W r (2796)-color anti-correlation reported by 
IZibetti et al] d2007l) . Both the W r (2796)-color correlation and 
the indistinguishable W, (2796 ) distributions in our sample are 
in conflict with the findings of Bor doloi et alj d201 ll) . who re- 
port eight times stronger W r (2796) within D = 50 kpc associ- 
ated with bluer galaxies. 

3. Including upper limits on W r (2796), we find that the 
significance level of the well-known anti-correlation between 
W r (2196) and impact parameter, D, is 8.2 a for this sample. 
The best parameterization is a log-linear relation, for which 
we find logW,.(2796) = (-0.019 ±0.002)£>+ (0.28 ±0.11). 
The scatter in this relation may be due to the redshifts, colors, 
or luminosities of the galaxies. Bifurcations in these galaxy 
properties result in at most a ~ 3.3 a correlation where higher 
^-band luminosity galaxies tend to have higher W r (2796) at 
a given impact parameter. Dividing the sample into quar- 
tiles based on the redshift, color, or luminosity of the galaxies 
and comparing the lowest and highest quartiles gives at best a 
~ 2.5 a correlation with B-band luminosity. 

4. The covering fraction profiles with projected distance 



from the galaxy, / Dm „ and f {D) , decrease with both increasing 
D and increasing minimum W, (2796) (see Figures [2] and [3]). 
High luminosity galaxies have higher covering fractions. In 
terms of f {D) , this difference is greatest at 100 < D < 200 kpc 
for absorbing CGM gas with W r (2796) > 0.1 A and at D < 
50 kpc for CGM gas with W r (2796) > 1 .0 A. There is no clear 
difference between the covering fraction profile of blue galax- 
ies versus red galaxies. The high redshift galaxies (z > (z)) 
have a higher covering fraction than the low redshift galaxies 
(z < {z}) at D > 50 kpc. Since our sample comprises a larger 
fraction of high luminosity galaxies at high redshift, this latter 
result may be tracing the luminosity dependence of the cov- 
ering fraction. A multivariate analysis would be useful for 
examining such possible cross-correlations in the sample. 

5. We determined the best-fit parameters B* and f3 for the 
luminosity scaling of the outer envelope for absorption, or the 
"halo absorption radius", R(L) = R*(L/L*)P, for both the B- 
and -bands as a function of W r (2796) threshold. For the 
B-band, the luminosity scaling remains roughly constant at 
f3 ~ 0.4 with increasing W cut . We find B* ~ 75 kpc for all 
but W cut = 0.1 A, which is on the order of ~ 90 kpc. The 
covering fraction inside R(Lb) decreases from f m = 0.89 for 
W C ut = 0.1 A to f m = 0.34 for W cut = 1.0 A. For the tf-band, 
the luminosity scaling flattens from j3 ~ 0.25 for W cut = 0.1 
and 0.3 A to f3 ~ 0.15 for W cut = 0.6 and 1.0 A, whereas B* 
decreases from ~ 75 kpc to ~ 60 kpc. The covering frac- 
tion behaves similarly to the B-band. Though the "outer en- 
velope" for absorption has a clear dependence on both lumi- 
nosity and W r (2796) threshold, we note that the scatter in the 
data in the W, (2796)-D plane, the behavior of the covering 
fraction profile f {D > with increasing D, and the relatively ex- 
tended upward uncertainties in B* all indicate that R(L) is not 
a well-defined quantity, but should be interpreted as a "fuzzy" 
boundary and that the "fuzziness" increases with decreasing 
equivalent width threshold. 

6. Dividing the galaxies by the median color (B-K = 1 .39), 
adopting W cut = 0.3 A, and examining the fit parameters to 
R(L) = R*(L/L*)P, we find that, in the B-band, blue galax- 
ies have a steeper luminosity scaling, j3 ~ 0.45, and larger 
covering fraction, f m ~ 0.9, than red galaxies, f3 ~ 0.3 and 
fmu ~ 0.75. The halo absorption radius for blue galaxies, 
B* ~ 90 kpc, is larger than for red galaxies, R* ~ 75 kpc. 
Our re sults are in conflict with the results of Bordoloi et al. 
J20T1 . 

7. We also examined the R(L) dependence in both the B- 
and -bands on galaxy redshift, bifurcating the sample at the 
median redshift z = 0.346, and adopting W cut = 0.3 A. In the B- 
band, we find B* ~ 90 kpc, j3 ~ 0.5, and f m ~ 0.8 in the low z 
subsample. The high z subsample has similar values forB* (~ 
85 kpc) and f m (~ 0.85), but a much shallower dependence 
on Lb/L b , with j3 ~ 0.3. In the /f-band, we find no discernable 
difference between the low and high z subsamples. The values 
for B„, f3, and f m for both subsamples are ~ 80 kpc, ~ 0.3, 
and ^0.8, respectively. 

8. We present the B- and /T-band luminosity func- 
tions, $(M), for subsamples split by W r (2796) < 0.3 A 
and W,(2796) > 0.3 A and by low z and high z- For 
all subsamples, the faint-end slope in the B-b and may be 
slightly flatter than the luminosity functions of Fab er et al] 
d2007h . $(M B ) for W r (2796) > 0.3 A at high z is sugges- 
tive of the turn down at Mb > ~20, which was reported by 
ISteidel. Dickinson. & Perssonl d!994l) . Comparing the low z 
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and high z subsamples, we found evidence for redshift evo- 
lution (3.6 er) in the fi-band for the galaxies associated with 
W r (2796) < 0.3 A absorption. In addition, for high z galax- 
ies associated with W,(2796) < 0.3 A, the mean color trends 
toward bluer galaxies, but at low z, the mean color trends to- 
ward redder galaxies. For W r (2796) > 0.3 A, the color se- 
quence along $(Mb) for the low redshift and high redshift 
subsamples shows no clear redshift trend. We also find a 6.3 a 
correlation between Mk and B-K such that lower luminos- 
ity galaxies in the AT-band tend to b e bluer (consistent with 
Stei del. Dickinson, & Perssonlll994T) . 

This work constitutes our first examination of the relation- 
ship of the Mg II absorbing CGM with galaxy properties us- 
ing MAGllCAT. In particular, we have examined the data 
by comparing direct observables and have avoided analyzing 
the quantities based upon converting or scaling the data to a 
preferred model of the CGM. Our aim has been to provide 
an unfiltered view of the Mg II absorbing CGM. Nonetheless, 
the simple bivariate testing we have conducted here does not 
provide full leverage for probing of the available data. In 
future work, we will apply multivariate analysis methods in 
which we also incorp orate mass estimates of the galaxies (see 
IChurchill et al.ll2012bl) . Mgll kinematics, and both low and 
high ionization absorption strengths of the CGM for these 
galaxies. We will also conduct comparative studies of the iso- 
lated galaxies and group galaxies in the catalog. 
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